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In the 1950's Hanbury Brown and Twiss showed that one could measure the angular sizes 
of astronomical radio sources and stars from correlations of signal intensities, rather than am- 
plitudes, in independent detectors. Their subsequent correlation experiments demonstrating 
quantum bunching of photons in incoherent light beams were seminal in the development 
of quantum optics. Since that time the technique of "intensity interferometry" has become 
a valuable probe of high energy nuclear and particle collisions, providing information on 
the space-time geometry of the collision. The effect is one of the few measurements in ele- 
mentary particle detection that depends on the wave mechanics of the produced particles. 
Here we discuss the basic physics of intensity interferometry, and its current applications in 
high energy nuclear physics, as well as recent applications in condensed matter and atomic 
physics. 
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I. INTRODUCTION 

> . 

, Hanbury Brown-Twiss (HBT) interferometry, the measurement of two identical particle correlations, has 
become a very important technique in particle and heavy-ion collisions, enabling one to probe the evolving 
geometry of the collision volume. Figures 1 and 2 illustrate the general idea of an HBT measurement: plotted 
in Fig. 1 is the two-particle correlation function, C(Qi nv ) - measured for pairs of 7r + as well as for pairs 
of 7T~ by NA44 for 200 GeV/A S on Pb at the CERN SPS - as a function of the invariant momentum 
difference <5i nv = [(pi — P2) 2 ] 1 ^ 2 of the two particles. The characteristic falloff distance Aq in momentum 

H— ' • 

I i of the correlation function is of order 50 MeV/c for pions; the length h/Aq, which is ~ 4 fm, is basically 
a measure of the size of the source of the final state pions, the size of the source when the pions no longer 

P ■ 

interact strongly with other particles. Also shown in Fig. 1, for comparison, is the correlation function 
for pairs of n + for 450 GeV protons on Pb, which, being broader, indicates a smaller source size. Figure 
2 similarly shows the correlation function for 7r + 7r + , ir~ir~, and K + K + pairs produced in collisions of Au 
on Au at 10.8 GeV/A measured by E877 at the AGS in Brookhaven, also as a function of the invariant 



momentum difference [||. 

In general, the two-correlation function is defined by 

c(q) = /<^;?*» (1) 

where np is the number of particles of momentum p measured in a single event, q = pi — $2, and the 
averages, denoted by (•••), are over an ensemble of events. The pairs in the numerator are taken from the 
same event, and the pairs in the denominator from different events. Usually, one also averages the numerator 
and denominator separately over a range of center-of-mass momenta P = p\ + P2 of the pair, an average 



* Lectures given at the XXXVII Zakopane School, June 1997. To be published in Acta Physica Polonica. 
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denoted here by {• ■ •}. As q becomes very large the correlations between the particles are lost, and the 
correlation function approaches unity. 

The basic issue I want to discuss in these lectures is how and why HBT interferometry works. The effect 
is in a unique class of experiments involving multiparticle correlations that are sensitive to the actual wave 
mechanics of particles as they stream out to the detectors. Normally, one imagines quantum mechanics as 
being important in high energy experiments only until the particles leave the interaction region; from then 
on one usually pictures them as little bullets on classical trajectories. (Quantum phenomena such as kaon 
regeneration and neutrino oscillations involve the internal degrees of freedom of the particles, rather than 
spatial amplitudes.) Considering the wave mechanics of the emitted particles in space and time is crucial to 
understanding questions such as why independent particle detectors give a greater signal when they are close 
together, corresponding to small q, than far apart. Further issues are: How accurately are the distances that 
are determined by the falloff of the correlation function related to the size of the system? In principal the 
correlation function at small momenta differences should rise up to 2 for a perfectly chaotic source. However, 
it only goes up to ~1.5 - 1.6 for the pion pairs shown in Figs. 1 and 2. What is the physics that reduces the 
correlation function at small momentum differences? What is the effect of final state Coulomb interactions 
on the measured correlations? 
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Fig. 1. Two-particle correlation function for ir + ir + and ir~ir~ pairs in 200 GeV/A collisions of S on Pb, and 
7r + 7r + pairs in collisions of 450 GeV p on Pb |Q. 
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Fig. 2. Two-particle correlation function for 7r + 7r + , 7r 7r , and K + K + pairs in collisions of Au on Au at 10.8 
GeV/A @. 

I will begin by describing the HBT effect in the simplest model of classical waves, and then discuss how 
one can understand HBT in terms of the quantum mechanics of the particles reaching detectors. Then I 
will turn to the nuclear physics applications and finally mention applications of HBT interferometry in both 
atomic and condensed matter physics. My aim here is to describe the physics underlying the HBT effect. 
For more detailed discussions of the current experimental situation in ultrarelativistic heavy-ion collisions 
and its theoretical interpretation, the reader is referred to, e.g., the reviews [3-8]. 



II. BASIC MODEL OF HBT INTENSITY INTERFEROMETRY 



HBT interferometry differs from ordinary amplitude interferometry in that it does not compare amplitudes 
(as in a Young's two-slit interferometer) but rather intensities at different points.]] The simplest picture of 
HBT interferometry, from which we can see the fundamental idea, is to consider two distant random point 
sources of light, a and b (of the same frequency), or more realistically for a star, a distribution of point 
sources, and imagine measuring the light falling in two independent telescopes 1 and 2; H see Fig. 3. The 
detectors are not connected by any wires. Assume that the sources are separated in space by R, the two 
detectors by d, and that the distance from the sources to the detectors, L, is much larger than R or d. 

Imagine that source a produces a spherical electromagnetic wave of amplitude ae lk ^ r ^ r "^ + ^ a /\f — f a \, and 
source b a spherical wave of amplitude /3e ifc ' r ~ ri> ' +l ^ i> /|r — where (f> a and 4>b are random phases (we ignore 
polarizations here) . Let us calculate the correlation of the electromagnetic intensities in 1 and 2 as a function 
of the separation of the two telescopes. The total amplitude at detector 1 is 

A x = \ ( ae ifcri " + ^ + ^e ifcrib+i0i> ) , (2) 

where r\ a is the distance from source a to detector 1, etc., and the total intensity in 1 is 



1 As we shall see below, there is a close connection between the two types of interferometry. 
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2 4- I «l 2 4- a * p e i (H r lb-ria) + 4>b-4>a) j_ ^,fl* p -j(fc(i-i[,-ri a )+0 i ,-0 a ) 



af3*e 



(3) 



with a similar result for I 2 . On averaging over the random phases the latter exponential terms average to 
zero, and we find the average intensities in the two detectors, 



1 



(Ji) = (h) = T2 ((M 2 ) + (|/?| 2 )) ■ 
The product of the averaged intensities (h)(l2) is independent of the separation of the detectors. 



(4) 




Fig. 3. Measurement of the separation of two sources, a and b, by correlation of intensities in detectors 1 
and 2. 



On the other hand, multiplication of the intensities I\h before averaging gives an extra non- vanishing 
term ~ (a* (3)(a(3*), and we find after averaging over the phases that 
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{hh) = (h)(h) + -^|a| 2 |/?| 2 cos(fe(ri Q - r 2a - r lb + r 2b )) 
|a| 4 + |/3| 4 ) + 2|a| 2 |/?| 2 (l + cos (fc(ri„ - r 2a - r lb + r 2b ))] . 



Then 
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cos {k(r la - r 2a - r u + r 2b )) 



(5) 



(G) 



For large separation between the sources and detectors (I> R), k(r\ a — r 2a — r± b + r 2b ) — > k(r a — fib) ■ 
(f 2 — f\) — R ■ (k 2 — fci), where fe, = kfi is the wavevector of the light seen in detector i. The correlated 
signal in Eq. (0) varies as a function of the detector separation d on a characteristic length scale 



d = A/<9, 



(7) 



where A is the wavelength of the light, and 9 = R/L is the angular size of the sources as seen from the 
detectors. Thus by varying the separation of the detectors one learns the apparent angle between the two 
sources, and with a knowledge of the individual wavevectors, the physical size of the source. 

If instead of two discrete sources, one has a distribution of sources, p(f), then averaging over the distri- 
bution, one finds that the correlation function measures the Fourier transform of the source distribution: 



C(d)-1 



d 3 rp{f)e l{ '' 1 ^ 2) - p 



(8) 



One important difference between astronomical observations and high energy physics is that the stars stay 
fixed, while in a collision, the system evolves a time scale of 10~ 23 to 10 -22 seconds, and thus one has to take 
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into account the changing geometry. As we will see, in high energy physics one measures not the Fourier 
transform of the distribution in space alone, but to good approximation the Fourier transform in both space 
and time. A second important difference is that in astronomy, in the absence of a knowledge of the distance 
to the source, one cannot measure the actual difference in direction of the wavevectors of the light in the 
two detectors, and thus one measures only the angular size of the source as seen from the detectors. In a 
high energy experiment, one can determine the wavevectors of the detected particles, and thus measure the 
absolute size of the source. 

To find an enhanced correlation at detector separation < X/8 it is not necessary for the two detectors to be 
wired together. One needs only to compare the data trains. Why, we may ask, do two independent nearby 
detectors produce extra signal? Essentially if the amplitude varies randomly then a positive fluctuation of 
the amplitude will produce a correlated increase in both measured signals, and vice versa for a negative 
fluctuation. For example, in black-body radiation, both the real and imaginary parts of the complex electric 
fields E ~ ae ik-r-iuit are Gaussianly distributed. For independent Gaussianly distributed real variables x 
and y, one finds simply that ({x 2 + y 2 ) 2 ) = 2(x 2 + y 2 ) 2 , so that 

(\E 1 \ 2 \E 1 \ 2 )=2({\E 1 \ 2 )) 2 , (9) 

while for a coherent source, e.g., a laser, (|i?i| 2 |ia| 2 ) ~ ((|i?i| 2 )) 2 . The extra factor of two is precisely the 
source of the HBT correlations, the enhancement that arises from the cosine term in Eqs. (||) and (||). 

Not apparent in the simple model above is how to deal with the time involved in making measurements. 
For example, how far apart can one shift the data trains in time in comparing the intensities in the two 
detectors and still find a correlation between the signals? I will return to these questions below. 



III. A BRIEF HISTORY OF THE HBT EFFECT 



The radar technology developed in the Second World War opened the field of radio astronomy in the 
postwar period, and soon led to the discovery of bright radio "stars" in the sky. One had no idea how big 
various sources, e.g., Cassiopeia A and Cygnus A, were, and the immediate problem was how to measure their 
sizes. The standard technique in use was Michaelson interferomctry, in which one compares the amplitudes 
of the light landing at two separated points, e.g., by converging the two signals using a lens and producing 
a diffraction pattern as a function of the separation of the points. From the structure of the diffraction 
pattern (on a distance scale X/9) one can determine the angular size of the source. Using this technique 
Michaelson measured the angular diameter of Jupiter's system of moons in 1891, and K. Schwarzschild first 
measured the angular diameter of binary stars in 1895. The resolution by amplitude interferometry at a 
given wavelength is limited by the size of the separations over which one can compare amplitudes. Were the 
radio sources to have had a large angular size, then one would only have needed a small separation of the 
two detectors. On the other hand, were the sources small, then it might have been necessary to separate the 
telescopes by distances too large, e.g., on opposite sides of the Atlantic, to be able correlate the amplitudes 
with the technology available in this period. This is the problem that the radio astronomer Robert Hanbury 
Brown at Jodrell Bank solved in 1949. His basic realization was that "if the radiation received at two places 
is mutually coherent, then the fluctuation in the intensity of the signals received at those two places is also 
correlated." PQ | Hanbury Brown then brought in Richard Twiss who had a more mathematical training to 
carry out the mathematical analysis of intensity correlations. 
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The first test of intensity interferometry was in 1950, when Hanbury Brown and Twiss measured the 
diameter of the sun using two radio telescopes operating at 2.4m wavelength (in the FM band) - quite a 
spectatacular demonstration of the technique. Their group then went on to measure the angular diameters 
of the Cas A and Cyg A radio sources, which turned out to be resolvable within a few kilometers. Since they 
could in fact have done Michaelson interferometry over such distances, Hanbury Brown described the intensity 
interferometry effort as "building a steam roller to crack a nut." [jllj Nowadays, Michaelson interferometry 
has completely replaced intensity interferometry in astronomy. In radio astronomy, amplitude interferometry 
is the basis of the Very Large Array (VLA) in Socorro, New Mexico, and the extended VLBI, in which one 
compares radio amplitudes in separated radio telescopes. The 10 m optical Michaelson interferometer on the 
Space Interferometry Mission satellite, to be flown in 2004, will be able to resolve objects to 5 microseconds 
of arc. 

Intensity interferometry actually has an intimate relation with Michaelson amplitude interferometry, as 
noted by Hanbury Brown and Twiss fl^| . Amplitude interferometry measures essentially the square of the 
sum of the amplitudes A\ and A 2 falling on detectors 1 and 2: 

\A X + A 2 \ 2 = \A,\ 2 + \A,\ 2 + {A\A 2 + A 1 A*). (10) 

The latter term in parentheses, called the "fringe visibility," V, is the interesting part of the signal. Averaged 
over random variation in the signal, V 2 is simply 

(V 2 ) = 2(\A 1 \ 2 \A 2 \ 2 ) + (AfAl) + {A\A* 2 ). (11) 

As one can see from the simple model above, Eq. (|J), the final terms vary rapidly on a scale of separations, 
d ~ A, the wavelength of the radiation, and average to zero. On the other hand, the first term in Eq. ( pi] ) 
is just twice the correlation of the intensities landing in the two detectors. Thus 

(V a )->2(hh); (12) 

the time-average of the square of the fringe visibility is proportional to the time-averaged correlation of the 
intensities. 

While it was well demonstrated both theoretically and experimentally that intensity interferometry worked 
for radio waves, which were commonly understood as classical fields, it was not obvious in the early 1950's 
that the effect should also work for light. Light being made of photons was more mysterious than radio 
signals made of classical electrical waves; the connections, now clear, were obscure at the time. Hanbury 
Brown and Twiss decided to test the idea for optics, with a simple tabletop experiment in which they used a 
beam from a mercury vapor lamp - a thermal source - and a half-silvered mirror to split the beam in two [ fujf . 
By measuring the intensity correlations between the two separated beams, they essentially compared the 
intensities at two different points in the unseparated beam, and by varying the relative path lengths between 
the mirror and the detectors they could vary the time separation, r, of the points. What they found was 
that while at large r there were no intensity correlations, the correlations increased with decreasing r. The 
characteristic timescale is the coherence time of the beam which, in this case is essentially h/T, where T is 
the temperature of the source. This experiment was the crucial demonstration of "photon bunching," i.e., 
that photons in a seemingly uncorrelated thermal beam tend to be detected in close-by pairs. Their results 
were greeted with great disbelief, and various experiments were done to disprove them. In the end Hanbury 
Brown and Twiss prevailed, aided by a particularly important paper by Purcell J14| which showed how to 
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understand the effect in terms of electric field fluctuations (see Eq. @) - and the field of quantum optics 
was born. 

Armed with the demonstration that intensity interferometry worked for light, Hanbury Brown and Twiss 
then went on to apply the technique to measure the angular size of the star Sirius (a Canus Majoris A) 
by studying optical intensity correlations between two telescopes (jl5). Since the telescopes required good 
light gathering ability but not great resolution, Hanbury Brown and Twiss were able to fashion a pair from 
five-foot diameter searchlights left over from the Second World War. The signals from the two telescopes 
were correlated electronically (although the actual physical connection is not needed to observe the effect). 
Figure 4 sketches their data for the correlation function C(d) — 1, divided by its value at d=0, measured as 
a function of the separation d of the two telescopes out to a maximum separation ~ 9 m. The data yielded 
an angular diameter of Sirius of 0.0068"± 0.0005" = 3.1 x 10~ 8 radians, a very impressive measurement of 
an object at a distance of 2.7 pc. The four data points shown were taken for a total of some 18 hours over 
a 5 month period, an indication of the poor viewing conditions. The dashed line is the expected curve for a 
uniformly illuminated disk of angular diameter, 0.0063". 



C(d)-1 



C(0)-1 



10m 



Fig. 4. Measurement of the angular diameter of Sirius pf 



This figure looks very similar in structure to the heavy-ion plots in Figs. 1 and 2. An important difference 
is that the actual HBT correlation seen here was just one part in 10 6 , a tiny signal above the background. 
What is the source of this difference? The question is whether all observed pairs of photons are "HBT- 
correlated"]^]; for example, if one takes a data train from 1956 in one telescope and compares it with a 
data train in the other from 1997 will one see interferometry? The answer is that the photons are in fact 
HBT-correlated only if they are emitted within a coherence time, the characteristic timescale in the original 
HBT tabletop experiment with an optical source. For a star the coherence time is T co h ~ 10~ 14 sec. On the 
other hand, the signal was studied over a band 5-45 MHz, corresponding to a binning time Tbi n ~ 10~ 8 
sec. Roughly the probability of observing an HBT-correlated pair of photons is ~ "Tcoh/Vbin ~ 10~ 6 . Figure 



2 Here we mean correlated in the sense that the photon pairs will produce an HBT effect at the detectors, measured 
by C, as opposed to the different question of the correlations in the beam produced in the source, as discussed in Sec. 
6. The language is potentially confusing, since photons that are maximally correlated at the source, e.g., in a laser 
beam, do not exhibit an HBT effect, while a thermal source, which is minimally correlated, produces the maximum 
HBT effect. 
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5 shows the region where photons produce an HBT signal in the plane of the times, t\ and t%, of detections 
in detectors 1 and 2. Below we discuss the analogous timescales in heavy-ion collisions. 




+ tbin 

Fig. 5. Region in the plane of the detection times t\ — where photon pairs produce an HBT signal. 

IV. QUANTUM MECHANICS OF HBT 

The simple derivation of intensity interferometry in Sec. 2 is entirely classical. How can one understand the 
effect from a quantum mechanical viewpoint? If we think of the sources a and b in Fig. 3 as emitting photons, 
we can identify four different processes, shown in Fig. 6: i) source a emits two photons, one detected in each 
detector, ii) source b emits the two photons instead; iii) source a emits a photon detected in 1 and b emits a 
photon detected in 2, and finally, iv) source a emits a photon detected in 2 and b emits a photon detected in 
1, the exchange of the previous process, iii. The first two processes are distinguishable, and do not produce 
any interferometry. They simply correspond to detection of the sources independently (the |a| 4 and |/3| 4 
terms in Eq. @). Only the latter two processes, iii and iv, which are quantum-mechanically coherent, give 
rise to interferometry. [Indeed, if we drop the terms proportional to |a| 4 and |/3| 4 , then Eq. (^) reduces to 
C{d ) = 1 + cos(fc(ri a — r-ia — ru, + ^2f)))-] Quantum mechanically, the HBT effect is a consequence of ordinary 
boson exchange, an effect included in the symmetry of the wave function of the pair of particles, e.g., for a 
pair of independent bosons in or thogonal states <p a and t)>p, 0(1,2) = (^(1)^(2) + ^(2)^(1)) /y/2. The 
effect is present for all pairs of identical bosons, including pions and kaons produced in high energy collisions. 

The detection of interferometry in particle collisions dates from 1962 when G. Goldhaber, S. Goldhaber, 
W. Y. Lee, and A. Pais [jl6| studied angular correlations of pions produced in pp collisions at the Bevatron. 
According to Pais [0, the idea of exploring interferometry in particle physics, although so similar to that in 
astronomical observations, was independently conceived. The method is now a standard technique in high 
energy collisions, from heavy ions [3-8], to meson-nucleon interactions |^8| , to electron-positron annihilation 
[p~9|| . As noted by Feynman [p(i| , the experiment done with electrons would yield intensity anti-correlations. 
However the effect is obscured by interactions among the fermions; electron pairs or proton pairs have a 
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repulsive Coulomb interaction which itself decreases the correlation function at small momentum differences 
(see Sec. 9), while neutrons at small relative momenta have significant final state strong interactions. 
Correlation studies of nucleon pairs produced in heavy- ion collisions are described in Refs. || and PH , and 
references therein. HBT interferometry is now being applied in study of boson atomic beams as well j22| , as 
we discuss in Sec. 10. 

Eventually, we will describe HBT measurements in terms of the two-particle correlation functions of the 
emitted identical particles. It is more intuitive, however, first to study the problem in terms of particle wave 
functions. To be specific we focus on pions, although the discussion is quite general. 




® ■ ID 

® > ID 

(iii) (iv) 

Fig. 6. The four independent photon emission and detection processes included in Eq. (|5|). 

How does one describe quantum mechanically the set of pions emitted emitted in a nuclear collision? Even 
in the best of all possible worlds - where one knows the exact wave function of the two colliding nuclei, and 
knows exactly how the quantum mechanical evolution operator does its job to produce the final system as 
a coherent superposition of well defined pure quantum mechanical states ^(l, 2, N) of N particles - the 
subset of pions emitted is described by a mixed quantum state. Quite generally, any subset of particles in a 
pure state is described by a mixed state, even, e.g., for the electron in the ground state of a freely moving 
hydrogen atom. The single particle density matrix for pions of a given charge at equal time is given by 

(7r t (r,t)7r(r , ,t))= / d 3 r 2 ■ ■ ■ d 3 r n ^*{r, r 2 , • • • , r N , t)V(r', r 2 , • • • , r N , t), (13) 

where tt(x) is the operator destroying a pion of the given charge at point x — (r,t). 

If 'J is a product of single particle wave functions, then (tp (x)tt(x')) factors into a product of single particle 
states, <jf{x)(j){x'). In general the single-pion correlation function does not factor, even when the pions are 
no longer interacting, but it can be represented as a sum over a collection of single particle states fa: 

{-k\x)-k{x')) =£>#(aO&(x'), (14) 

i 

where the fi give the probability of finding the pion in single particle state i. For example, the probability of 
finding a pion at point 1 is given by ^ fi\fa(l)\ 2 . Only if the pion part of the state \& factors out in the form 
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of a product of the same single particle states for all the pions - a Bose-Einstein condensate - will single pions 
be in a pure state. The mixed single pion state always has finite entropy, — (/* In /, — (1 + fi) ln(l + /j)). 

The closest one can come to describing pions as little bullets is to picture the single particle states making 
up the mixed ensemble of pions as wavepackets with almost well defined momenta and energies, limited by 
the uncertainty principle. We picture the collision volume as made up of many sources of pions; whether 
the sources are fragmentation of strings, or in the language of low energy nuclear physics individual nucleon- 
nucleon collisions, the sources are localized to within a distance R which is less than the size of the entire 
collision volume, and the emission process is temporally localized to within a time r. Thus the individual 
components of momentum, p, and energy, e p , of the emitted particles are uncertain to within 

Ap a >h/R, a = x,y,z 

Ae p >h/r. (15) 

A pion nominally of momentum p emitted from a source at the origin in space and time would have an 
amplitude to have four-momentum q — (e q , q) that is roughly Gaussian, 

A(q) ~ e -(9-P) 2 fl 2 /2 e -(e 9 -e P ) 2 r 2 /2 5 (16) 

and in space-time the particle would be described by a wavepacket 

d 3 q e iqx 
(2^)3 2T q 



M*)= I 7^3— A (l)- (17) 



How does the packet evolve in time after leaving the source? Assume that the collision takes place, and 
the particles emerge, into vacuum. (In real life secondary scattering on other atoms in the target and also 
scattering in air are both very important effects, to which we return in Sec. 8.) The transverse spatial 
spread, perpendicular to p, is determined by the uncertainty in transverse velocity, which for for relativistic 
particles (s p s=s p) is: 



Av± _ = ^ ~ -L. (18) 
e p pR 



The spread in longitudinal velocity is 



Av l = a( P -±)=^~' (19) 
\s p J yep tpR 

where 7 = e p /m is the Lorentz factor, and m is the pion mass. The presence of the factor j 2 reflects the 
fact that the more relativistic the particles become the closer to the speed of light the packet moves, with 
vanishing spread in longitudinal velocity. The transverse size of the wavepacket after travelling a distance L 
is thus <~ L/pR, while its thickness is ~ L/j 2 pR. 

To be specific, consider a 1 GeV pion (7 = 7) produced within an initial radius of 10 fm travelling to a 
detector 10 meters from the source. Then the transverse spread of the wavepacket is about 20 cm, while the 
thickness of the packet grows to half a centimeter. Pions emerge in rather extended pancake-like states. The 
characteristic time for such a pion wavepacket to cross a point at the distance of the detector is <~ 10~ n 
sec. Smaller source sizes lead to even larger spreads in the pion wavepackets. Photons would have a similar 
transverse spread; however, the thickness of the wavepacket would not grow. 
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V. DETECTOR RESPONSE 



Let us now turn to the question of the detection of particles in such wavepacket states by a magnetic 
spectrometer. In the "bullet" picture, the particle travels along a classical trajectory, excites an atom in the 
detector, which determines the direction of its momentum, and begins to make a track (or in emulsion makes 
a spot on photographic film); the particle continues on, producing a curved track in the spectrometer, from 
which one deduces the magnitude of its momentum. But, in fact, the particle has a distribution of probability 
amplitudes, wherever its wavepacket <fi is non-zero, of where it makes the first spot in the spectrometer. A 
particle of mean momentum p can, because of the transverse momentum uncertainty, be detected anywhere 
within the range of momenta Ap± around p contained in its wavepacket. After the initial collision in the 
detector, a much more narrowly focused wave packet emerges. (If one does not actually do the measurement, 
the state that emerges is a mixed state corresponding to all possible points where the incident packet can 
excite a detector atom.) The narrowed packet continues through the magnetic field, and the subsequent 
collisions it makes selects its energy. One measures momentum by measuring a sequence of positions in the 
detection. 

Consider the measurement of the momentum of an incident particle in state <t>i(x). The probability that 
it makes the first single excitation of an atom in the detector at point A and continues to produce a track 
corresponding to measuring its momentum to be k is then 



Here the detector atom response function, s, is localized in r and r' about the atom at A; it also extends 
over an interval in t — t' which is the characteristic time that the detector analyzes the amplitude and 
phase variation of the incident wave, of order 10~ 16 sec ~ h/(10 eV), the inverse of a characteristic atomic 
excitation energy. This time is much shorter than the typical nanosecond scale resolution time of a detector, 
the time it takes to build up a million-fold cascade of electrons. 3 More generally, the probability of detection 
of a pion of given charge at A is given by 



where (tp (x)tt(x')} is the single-pion correlation function. 

Consider next the detection of two independent particles. The question is why, if one detector lights up, do 
nearby detectors tend to have greater probability to light up than detectors further away - the HBT effect. 
Suppose first that the two particles are incident in orthogonal wavepackets <f>i and <pj. The symmetrized 
two-particle wave function is (j>(r, f \t) = ((t> i (r f ,t)(j) j (r',t) + (t>i(r f ',t)<l> j (r,t))/V2- The probability of a joint 
detection by a detector atom at A of a pion that continues to make a track corresponding to momentum fc, 
and by a detector atom at B of a second pion that continues to make a track corresponding to momentum 
k ', is then given by 




(20) 




(21) 




3 A detailed discussion of the role of the detectors is given in J. Popp's thesis |23|], and in Ref. 
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= P~ k (A; i)P % , (B- j) + P- k (A; j)Pp{B- i) 

+ A k (A; i,j)A k ,{B- j, i) + A % (A- j, i)A k , (B; i, j), (22) 

where 

At(C;i,j)= [ dxdx"e lk ^- x "h c (x,x")(j)*{x)(t) J {x"). (23) 



The first term in Eq. (|22j) is the probability of the particle in state i being detected at A times the probability 
of the particle in state j being detected at B, and the second is the same with i and j interchanged. These 
are the normal terms. 

The final two terms in Eq. (^2|) are the enhancement of the detection probability - the HBT effect. As 
we see, in order to have enhancement, it is necessary that both wavepackets, i and j, overlap in the detector 
at A during the time that the detector is doing quantum mechanics on the incoming system, and similarly 
that they also must overlap in the detector at B (but note that the wavepackets do not have to be present 
in both detectors simultaneously) . The presence of the wavepackets simultaneously in each of the detectors 
is the reason the detectors produce more signal when they are close to each other .0 

The maximum transverse separation d of the detectors that will produce an HBT signal is essentially the 
transverse size of a given wavepacket, L/pR, the wavelength divided by the angular size of the source as seen 
from the detectors, where again L is the distance from the source to the detector, R is the size of the source, 
and p is the average particle momentum. The scale of relative momenta q for which one finds a signal is 
q/p = d/L ~ 1/pR, and thus q ~ 1/-R, the standard HBT result. As we see from this argument, the HBT 
effect directly measures the width of the wavepackets at the first detection. This width is in turn determined 
by the uncertainties in the momentum distribution at the time that the wavepacket is no longer affected by 
strong interactions with the other particles in the collision.. 



Generally the amplitudes, Eqs. ( |16| ) and (17), of the wavepackets, A(q), vary slowly in q on an atomic 
scale, ~ 1 KeV/c. In this case the overlap integral ( p3| ) is not sensitive to the detection time scales, ~ 10~ 16 
sec, and one finds a correlation function: 

C(v v')-l+ \EiMM^M)l m) 

%P)_1+ WFE^W' (24) 

where fi is the single pion probability in the ensemble, Eq. (|l4|). 

However, if a particle is delayed in emission by more than the detection time scale, the HBT correlations 
between that particle and one produced directly will be suppressed. A simple example in heavy-ion collisions 
of this effect is in the correlation of n~ produced in A decay, A — ► n~ + p, with directly produced pions. 
Because the A travels more slowly than a directly produced pion of the same rapidity as one emitted in 
the decay, the pion from decay will lag the directly produced one by a time At. To estimate this effect, 
we note that a ir~ emitted in the forward direction has rapidity y° ~ 0.67 in the A frame, and that a 
A of rapidity y travels on average a distance T\smhy before decaying, where ta is the A lifetime. Thus 



4 Imagine that the detection at A occurs before the detection at B. Then one may ask how both original wavepackets 
can be at the second detector, since the first detection "reduces" the wavepacket of the detected particle. From this 
point of view the amplitude for detection of a particle at B is proportional to the amplitude, <pj(B), for it to be in 
state j at B times the amplitude, cj>i(A), for the other particle in state i to have been detected earlier at A, plus the 
same product with i <-> j. The resulting joint probability is the same as Eq. (p^). I thank J. Walcher for raising this 
question. 
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At = TA/(coshy + sinhy/ tanhyJJ), which for a A of typical rapidity 3 is ~ 0.037ta = 9.7 x 10~ 12 sec, much 
longer than the detector timescale. Pions emitted in other than the forward direction will have an even 
greater time lag. 
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Fig. 7. a) Two stars, S, and T, along nearby lines of sight from the earth; b) schematic of HBT measurement 
of correlated intensity from the two stars; c) schematic of HBT measurement a bright star surrounded by a 
halo of dim stars. 




RHIC 




(cern) 

Fig. 8. Detection of the HBT effect between a pion from RHIC and a pion from CERN. 

As the above discussion makes apparent, the enhanced signal is not a consequence of special preparation at 
the source, such as a particle of a given momentum inducing emission of other particles of similar momentum, 
as in a laser. (In fact, a coherent source such as a laser would not give an HBT signal.) Clearly, in the stellar 
case there can be no such connection between emission processes on opposite sides of the star, and yet photons 
from opposite sides give an HBT enhancement. The effect is a property of the detection. Furthermore, if two 
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stars, S and T, which are at very different distances from the Earth, are approximately along the same line 
of sight, Fig. 7a, one will see HBT correlations between a photon from S and a photon from T. Of course 
the emission times have to be different in order that the two photons arrive at each of the detectors at the 
same time. (Even in measuring a single star one correlates earlier emitted photons from the stellar rim with 
later emitted photons from the front surface.) The HBT correlations between photons from the two sources 
as a function of detector separation, d, sketched in Fig. 7b, would have an oscillatory term characteristic 
of the angular separation of the two stars (cf. Eq. (|^), modulated by a broader Gaussian characteristic of 
the angular size of an individual star. By comparison, a single bright star, T, surrounded by a halo of dim 
stars (an average over a distribution of stars S), would yield a correlation function, Fig. 7c, with a rise 
proportional, for to the inverse of the angular size of the central star X/d,2, plus a much more narrow 

rise, for d^d\, inversely proportional to the angular size of the halo. As an exercise, sketch the correlation 
function produced by two lasers which are mutually incoherent, rather than the two stars. 

That the HBT effect is a not a property of production, but rather is a property of detection is well 
illustrated by the following experiment. Imagine, as in Fig. 8, that RH1C, the BNL heavy-ion collider, sends 
a pion into a far away particle detector array, and that by using the pionic analog of a half-silvered mirror, 
we reflect a pion from CERN into the same detector array, so that it arrives at the same time as the pion 
from RHIC. As long as the wavepackets of the pions overlap sufficiently that the interference amplitudes A 
at the detector atoms are non-zero, one will see an HBT enhancement! As this experiment illustrates, the 
HBT effect does not depend on the history of the particles. What matters is the form of the wavepackets 
when they arrive at the detectors. 



VI. CORRELATION FUNCTIONS 



Let us now turn to describing the HBT effect in terms of the correlations of the sources of pions in a 
collision. To generalize Eq. (|22|), we may write the joint probability of detection at the detectors at A and 
B in terms of the two-pion correlation function {t:* {x)tt {x")n{x"')n{x')) as 

p%,p (AS) 

dxdx'^ x - x 'h A {x,x') / dx"dx"'e ik '( x "- x "'h B (x",x"') 

x (ii\x)^{x")it{x"')it{x')). (25) 

The two-pion correlation function (7r' (1)7T (2)7r(3)7r(4)), where the integers stand for space-time points, is 
the amplitude for starting in a state of the system, removing a pion at 4, then removing a second pion at 3, 
adding a pion back at 2, adding another back at 1, and returning to the initial state. 

Consider first a system of N free bosons that is completely Bose-Einstein condensed, or the photons 
produced by a laser. For such a system all the particles, or photons, are in the same particle state 4>{x). 
Then the correlation function completely factors. The amplitude for removing a particle at a given point 
4 is simply proportional to the wave function cj> at the point, while the amplitude for adding a particle is 
proportional to <fi* at the point. In this case the single-particle correlation function is 

(7rt(l)7r(2))=^* (1)^(2) (26) 

and 
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(7rt(l)7rt(2)7r(3)7r(4)) = A^* (1)0* (2)0(3)^(4) 
= < 7 rt(l) 7 r(4))( 7 rt(2) 7 r(3)). 



(27) 



Such a source is co/iereni.[j 

The correlations of particles from a thermal source, e.g., a black-body cavity, are quite different. For free 
particles, the mean number of particles of energy e is given by (n e ) = 1/ (e^( e— ^ — 1), while the fluctuations 
are given by 

<^)-(n £ ) 2 = <n £ )(l±(n £ )), (28) 

where the upper sign is for bosons and the lower for fermions. (The fermion result follows trivially since 
n can only equal or 1.) Translated back into a statement about the correlation function, one finds the 
expected factorized form for a thermal ensemble: 

(^(1)^(2)^(3)^(4)) = (^(1)^(4)) (^(2)^(3)) ± (^(1)^(3)) (^(2)^(4)), (29) 

where 7r here represents either a Bose or Fermi field. This equation says that if one removes two particles at 
4 and 3, one can come back to the same state by either replacing the first with a particle at 1 and the second 
at 2 (first term), or the first at 2 and the second at 1 (second term). In the boson case, when 3=4 and 1=2, 
the amplitude for removing two particles is just 2(7r'l'(l)7r(4)) 2 . The extra fluctuations - the factor of 2 here, 
or more generally the second term on the right side of Eq. ( p9j ) - are the source of the HBT interferometry 
effect. The maximum HBT effect occurs when the correlation function factorizes in this fashion. Then one 
describes the source as chaotic. 

The source need not be thermal to factor this way. The result ( |2^ ) always holds for non-interacting 
fermions, while for bosons it is sufficient that no single particle mode i be macroscopically occupied, i.e., 
that all (m) are <C 1. The basic reason one expects the correlation function in heavy- ion collisions to factorize 
as in Eq. ( |29| ) is that the pions undergo considerable rescattering in the hot environment of the collision 
volume; the key is the destruction by rescattering of phase correlations among the pions from the production 
process.^] 

It is useful to relate the pion correlation functions to correlations of the sources of the pion field. The 
freely propagating field ir(x) measured at the detector is produced according to the field equation 

{U 2 -ml)n{x) = J(x), (30) 

where J(x) is the source of the field at the last collision. Then 

{ ( ii 3 k p ikx ( i 

n(x) = / dx'D(x,x')J(x') = / — / dx'e~ lkx J(x'), (31) 

J J (27r)' 3 2ie k J 

where D is the free pion Green's function, and the latter form holds in the far field. The single pion correlation 
function is related to the source correlation function by 



5 Taking Poisson statistics for the distribution of the number n of photons in a laser beam more carefully into account 
leads to the same result @, since (n(n — 1)) = (n) . 

6 A simple example that violates Eq. (^9|) are the correlations among pions radiated by a weakly interacting gas of 
nucleons. 
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k\x)h{x'))= / dx"dx"'D*(x,x")D(x',x"')(j\x")J(x"')). (32) 



Generally (tt^ (x)tt(x')) can be written in a factorized form as a sum of wavepackets as in Eq. (|l4|). 
The measured singles distribution of pions is then given in terms of the source correlation function by 

d 3 n d 3 n 



Bp- 



d 3 p d 2 p±dy 2 
and the measured pair distribution by 
, d 6 n 1 



dxdx'e' p{x - x ' ) (j\x)J(x')), (33) 



Zp d 3 pd 3 p' 4 



dxdx l e lp{x - x ' ) dx"dx l "e lp ' (x "- x "'\j\x)J^{x")J{x"')J{x 1 )). 



(34) 



In the following let us assume a chaotic source, so that 



{J\x)jHx")J{x"')J{x')) 
(j\x)J(x'))(j\x")J{x'")) + {J\x)J(x"')){.P{x")J{x'))- (35) 



then the pair distribution function becomes 
_, d e 



£p£p d 3 pd 3 p 



- = \ f dxdx l dx"dx'"e l ^ x - x ' ) e lp(x ''- x ''' ) {j\x)J{x l )){.P{x")J{x'")) 

e ip(x-x') e ip'(x"-x"') , e ip(x-x"') e ip' (x" -x') 



(36) 



We see that the HBT correlation function, defined by 



d*n/d 3 pd 3 p> 

W {d 3 n/d 3 p){d 3 n/d 3 p'Y K ' 

where q = p — p' (and with implied separate averages over the center-of-mass coordinates of the numerator 
and denominator, cf. Eq. (|l|)), measures the structure of the current-current correlation function. Note that 
the information it provides is on the nature of the source of particles after the last strong interactions, when 
the particles begin to stream freely towards the detectors. 

The correlation function ( J'(x) J{x')) contains two length, and time, scales. The range of the center-of- 
mass variables, X = (r + r')/2, (t + t')/2, are on the order of the size of the collision volume, R ~ 10 fm, 
and the duration of the collision, r, also on the order of 5 - 10 fm/c. On the other hand, the dependences in 
x — x' measure the space-time extent of the region in which the phase at a point x' is coherent with the phase 
of the current at x, a region of size, £ c and r c in space and time. The lengths £ c and t c , which determine 
the falloff of the singles distribution, Eq. (^), are typically on the order of one fm in space and one fm/c in 
time, much shorter than the range in the center-of-mass variables. 

Such behavior is illustrated by the factorized form for (, (x)J(x')} , 

(JHx)J(x')) = e -(-+-') 2 /8« 2 e -(*+*') 2 ]/8r 2 5(a . _ ^ (38) 

The singles distribution, from Eq. (|33|), is then 

d 3 n f 

Sp tfp~~ J dxeWX ^ ( 39 ) 
and the two-particle correlation function, assuming a chaotic source, is 
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C(a) = l + e-« 2fl2 e-«° V {d 3 n/d 3 Kf 

W (d 3 n/d 3 (K + q/2))(d 3 n/d 3 (K - g/2))' 1 j 

where K = (p + p')/2. As we see from this equation, the length measured in HBT is modified from the 
length, R , governing the center-of-mass behavior of the current-current correlation function, (J^ (x)J(x')) , 
due to the final factor, the ratios of the singles distributions. 

A particularly simple and illustrative model is the following. Assume that the particle production is 
described by a distribution of sources of size R s , r s , at space-time points x s , each producing pions in 
wavepackets of mean momentum p , 

(41) 

with probability f{p)- The spread in momenta in the individual states 4> P (x) is of order h/R s . If the 
sources are Gaussianly distributed in space and time over a region of size Ro,tq, then from Eq. (|l4j) the pion 
correlation function is, 

(^(xHx')) ~ J 7^3 /(P) J d^e-^h-^o^x-x^x' - Xs ). (42) 

Carrying out the x s integrals explicitly, one readily finds that this model is the same as that with a source 
of the form (|38|), where 



R' = R£+ Rt/2, r z = rt + ii/2, (43) 



and 



9(*) = J ^f(P>- WX ^ x2/4R2 °- (44) 

The details of the individual wavepackets are all subsumed in the current-current correlation function. As 
this model illustrates, the length scale describing the ccntcr-of-mass of the current correlation function is the 
size of the distribution of sources, plus a correction from the size of the individual sources. 

In fact, this latter correction is countered by the correction from the singles distributions in Eq. (flO]). If 
we assume, solely as a mathematically simple example, that d 3 n/d 3 K ~ e^^ K ^ , then 

c( q ) = i + e -^^+(Ri-e)m e - q °\ri + r-,2)^ (45) 

The net deviations of the measured scales R and r from the scales of the source distributions Rq and tq are 
of order of a few percent at most. 

More generally one can define a total-momentum dependent pair source function, 

S P (X) = / dxe- iPx ' 2 {j\X + x/2)J(X - x/2)), (46) 



where P = p + p' is the total four-momentum of the pair. For the simple example (|38|), we have Sp(X) — 
e -x /2R g(p/2); note the relation to the description ([l2|) in terms of wavepackets produced by the source. 
In terms of S, the HBT correlation function becomes [p5| , 

C(q ) = l + \SdXe^S P (X)? 

° w + jdxs P+q/2 (x)fdxs P _ q/2 (xy 1 > 
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This equation relates the HBT correlation function to four-dimensional Fourier transforms of the source 
function Sp(x). Compare with the result (|24|), which gives the correlation function in terms of the Fourier 
transforms of the wavepackets making up the pion distribution. 

The simplest approximation is to take £ c = t c = 0, or equivalently, to ignore the P dependence in Sp(X). 
Then the current-current correlation becomes a function of only one variable: 



{j\x)J(x')) S(x)S(x-x') 



(48) 



This approximation is excellent for stars, where the emission of a photon is coherent on the order of an 
atomic scale, while the size of the star is ~ 10 11 cm. Neglecting the correlation length is not as good an 
approximation in a heavy-ion collision. 

With the neglect of the finite size of the correlation lengths, we find from Eq. (|33|) that 



(49) 



i.e., the singles distribution is flat in momentum. Furthermore the pair distribution (assuming a chaotic 
source) becomes 



£p£p d 3 pd 3 p> 



If 2 1 f 

— / dxS(x) + — / dxS(x)e 



i(p—p )x 



and 



C(q) = 1 + 



|J dxS(x)e lqx 
\J dxS(x)\ 2 



(50) 



(51) 



VII. PARAMETERIZATIONS OF DATA 

The most straightforward way to analyze HBT data is to parametrize the correlation function C(q) as a 
Gaussian in q. Expanding C(q) in Eq. ( |47| ) for small q to second order we find 

C(q) = 2 - q^dX^Xv) - {X„){X V )) 

+ <F<f({xpX v )-{x li ){x v )) + >--, (52) 

where 

_ f dxdX9(x,X)(jHX + x/2)J(X - x/2)) 
W**" JdxdX{Ji(X + x/2)J(X-x/2)) ■ 1 ' 

The terms in x in (|52|), which come from expansion of the denominator in Eq. (^7j), are of relative order 
(£ C /R) 2 , (t c /t) 2 , and are often neglected in the interpretation of the data, although, as mentioned, they can 
modify the extracted sizes by a few percent. Dropping these latter terms we have 

C(q) = 2 - q»q v {(x^X v ) - (x M )(x v )) +■■■ (54) 

where here and below (0(x)) = f dxSp(x)9(x) / J dx' Sp(x'). This form suggests a parametrization p6| , p7f , 

C(q) = 1 + Xe-^i^^^-^^, (55) 
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where we also introduce the chaoticity parameter A; for a completely chaotic source the correlation function 
rises up to 2 as q — ► 0, and thus A = 1, while for a completely coherent source, such as a laser, A = 0. 

Various increasingly sophisticated versions of this parametrization have been adopted. The simplest is to 
write 

C(q) = l + Xe-Q^ R \ (56) 

where Qi nv is the invariant momentum difference of the two particles, Qf nv = (pi —p-i) 2 — {spi ~£p 2 ) 2 - Results 
of such a single-size analysis by NA44 for pairs of ir + and pairs of tt~ produced in 200 GeV/A collisions of S 
on Pb at the SPS § are shown in Fig. 1, and by E877 for pairs of tt + , pairs of 7r , and pairs of K + produced 
in collisions of 10.8 GeV/A Au on Au at the AGS Q in Fig. 2. This parametrization corresponds to the 
assumption that {x^x v ) — (x fl )(x u ) — g^ v R 2 . Since the sign of the contribution of the time-time component 
should be the same as the space-space components a somewhat better single-size parametrization would be 
to assume that {x^x v ) — (x^)(x w ) = S^R 2 . Then 

C(q) = l + \e-^ 2 +i° 2 ) R2 ; (57) 

cf. Eq. ©. 

A second level of approximation is to distinguish the space and time dependence of the evolving system, 
taking a spherical fireball in space, so that 

C(q) = l + Xe-^ 2R2+ ^ T2 l (58) 

The time r is essentially the duration of the collision: t 2 = (t 2 ) — (t) 2 , and R the radius of the collision 
volume: R 2 = (r 2 ) — (r) 2 . 

The next level is to try to take the evolving geometry into account, including non-sphericity of the 
source and possible flow effects. Consider a pair of particles of total three-momentum P and relative three- 
momentum q. Since q ■ P = (p — p') ■ (p + p') = p 2 — p' 2 = 0, we have 

q° = q-P/P° = q- v, (59) 

where v is the velocity of the center-of-mass of the pair of particles, P/P . Then q^x^ = q ■ (r — vt), and 

qt><n{xnx v ) - (x^)(x„)) = ((q ■ (r - vt)) 2 ) - (q ■ (r - vt)) 2 . (60) 

Let us erect a three dimensional coordinate system in which the longitudinal direction is along the beam 
axis, the outwards axis (the x direction) is along the transverse component of P, and the third, or side, axis 
is in the y direction. In this frame v y vanishes. (Note that this coordinate system varies with the pair of 
particles studied.) The ensemble of events is symmetric under y — > —y, so that cross terms involving a single 
y vanish. However, the cross terms (zt) — (z)(t) and (xt) — (x)(t) are generally non-zero, and we find 

= 9o 2 ut (((x - v x t) 2 ) -{x- v x t) 2 ) + q 2 idc ((y 2 ) - (y) 2 ) 

+ 2g out giong(((a; ^ v x t)(z - v z t)) - (x - v x t)(z - v z t)), (61) 
a form which suggests a parametrization of the correlation function in terms of four radii pa] , 
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where the parameters have the interpretation 



(62) 



R 



2 

side 



R 



long 



= ((a; - v x t) 2 ) - (x - v x t) 2 , 

= (y 2 ) - (y) 2 , 

= ((z-v z t) 2 ) - (z-v z t) 2 , 



R l\ = (( x ~ v x t)(z - v z t)) -(x- v x t)(z - v z t). 



(63) 



Note that R%, although written as a square, need not be positive. It is often convenient to analyze data, 
pair-by-pair, in the "longitudinal center-of-mass" frame, in which P z = 0; then q out = q /v, and i?i 2 ong 
reduces to (z 2 ) ~ (z) 2 . 



The three-dimensional parametrization, Eq. (63), is commonly used in interpreting present HBT mesure- 
ments. Typical three-dimensional analyses of correlations of pion pairs are shown in Figs. 9 and 10, in Fig. 9 
correlations of of tt + it + from 200 GeV/A S+Pb collisions studied by NA44 ^9|, and in Fig. 10 correlations 
of 7T+7T+ and tt-tt- from 10.8 GeV/A Au+Au collisions studied by E877 Q. 

Considerable information on the development of the collision volume, e.g., flow |n|, can be extracted 
from the three-dimensional analyses. The experimental dependence of the two-particle correlations on the 
momenta of the particles indeed indicates that the systems are expanding. For detailed discussions of the 
physics extracted from recent experiments, see, e.g., Refs. (p2[ and p3[. 
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Fig. 9. Three-dimensional fittings of tt + tt + correlation functions, from NA44 [ p9| , the upper left panel as a 
function of g s ide, the upper right as a function of gout, and the lower panel as a function of gi on g- 
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Fig. 10. Three-dimensional fittings of ir + ir + and ir tt correlation functions, from E877 
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Fig. 11. Two-dimensional fittings of K + K + (left panels) and K K (right panels) correlation functions, as 
functions of gtransv and giong ■ From Ref . [pf . 
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Figure 11 shows NA44 data on K + K + and K~K~ correlations, projected as functions of qi on g and the 
component qtransv of q perpendicular to the beam axis || . Note that the chaoticity parameter A is somewhat 
larger for kaons than pions, a point which examine in the following section. 

VIII. SOURCES OF CHAOTICITY A < 1 

The chaoticity parameter, A, is generally found experimentally to be less than one, a reflection of intrinsic 
physical effects as well as experimental difficulties. The most fundamental effect would be that the source 
exhibits a level of coherence, the situation in a laser, or a form of pion or other boson condensate. HBT 
measurements of pions produced from a disordered chiral condensate in an ultrarelativistic heavy-ion collision 
would also show a reduced A jj4|. However, as we noted above, rescattering by other particles in the collision 
volume tends to destroy phase correlations from the production process. Another example is the MIT atom 
laser ]3o|l where magnetically trapped and evaporatively cooled sodium atoms are extracted in coherent states 
from a Bose-Einstein condensed system; since the extracted atoms do not exhibit an HBT effect, A would 
be zero. (See the discussion of HBT in atomic beams below.) 

Even if the source is completely chaotic, measurements do not necessarily give A =1. The first reason is 
the simple but important problem of contamination of the sample from misidentification of particles, e.g., 
an e~ or K~ as a 7r~, so that one includes pairs of non- identical particles in the data set. 

A second stems from unravelling the effects of Coulomb final state interactions between a pair of identical 
charged particles. The point is that the methods of removing effects of Coulomb interactions, which we 
discuss in some detail in the following section, become more uncertain the smaller the relative momentum 
difference, leading to uncertainty in A. 

The next physical effect reducing A is the production of pions from long-lived resonances. Such pions 
appear to come from sources of large radii, which would give an HBT enhancement only at very small q. 
Indeed, perhaps half of the pions produced in an ultrarelativistic heavy-ion collision come from resonances, 
rather than being produced directly. Pions from short-lived resonances, e.g., from p — > irir or from A — * Nir, 
are produced well within the collision volume and are not an issue. On the other hand, the long-lived 
resonance r/ has a lifetime of order 1.2 A/c, and the 3-7T into which it decays would appear to be produced 
at a relatively enormous distance of order A from the collision volume. The to goes some 24 fm, the rf some 
800 fm, etc. The result is that the collision volume is surrounded by a halo of pions from resonances. 

A small chaotic source would lead to a broad HBT correlation function, while a very large source would 
lead to a correlation function with a sharp rise only at small relative momenta. Now if one has both a small 
and a large source, e.g., a partially transparent cloud in front of the sun, one sees a combination of both, 
as illustrated in Fig. 7c, where the width of the bump closest to the origin is inversely proportional to the 
size of the large source (the cloud), and the width of the broader bump reflects the size of the small source 
(the sun). Figure 12 illustrates how pions from different resonances contribute to the correlation function, 
here as a function of the out component of the momentum difference, in a central CERN S-Pb collision [^G| . 
The estimate is based on an RQMD simulation of the collision. Note the rise at very small relative momenta 
from long-lived resonances. f\ However, unless one is capable of resolving the little peak at small q one would 



7 The calculation assumes that the pions from resonances are described by the same factorized form of the two-pion 
correlation function as the pions emerging directly from the collision volume. However, the pions from longer-lived 
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deduce that the data goes to a value of A less than one. The effect on the observed A is a reduction of ~ 
30% for pions, and ~ 10% for kaons, since a smaller fraction of kaons are produced by long-lived resonances. 
Since the production of resonances falls off with increasing transverse momentum, one finds in fact that the 
contribution of pions from resonances to the HBT signal decreases as p± of the pions studied increases. 

In the simple model of HBT described above we assumed well defined wave functions propagating through 
vacuum. But in reality the particles propagate through of order 1 mm of target and then, e.g., in NA44, 
through 15 meters of air. Let us consider the effects of secondary scattering in the target and the intervening 
air. From a quantum mechanical point of view scattering from the target and air atoms changes an initial 
pure state wave function of a particle into a mixed quantum state: when a pure wave function hits the atoms 
in the target or air, it generates, a la Huygens, a beam of secondary wave functions; because the atoms are 
disturbed, the secondary waves become incoherent with the initial wave. The secondary interactions produce 
many small angle scatterings of the particles, which have the observational effect of distorting the correlation 
function jj7|. 
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Fig. 12. Estimated contributions of pions from resonances in HBT compared with NA44 data 0. From fief . 

To see the effects of small angle scattering, it is adequate to describe the interaction of a high energy particle 
with an atom of charge Z by a screened Coulomb potential, V(r) — Ze 2 e~ r / a /r, where a — 0.8853CLQ/Z 1 / 3 
is the Fermi-Thomas radius of the atom and clq is the Bohr radius, and neglect elastic scatterings due to 
strong interaction with the nucleus. The differential cross section of an incident particle of momentum p is 
given in the Born approximation by 

da_ = Ap 2 Z(Z + l)a 2 

dQ {q 2 + l/a 2 ) 2 ' 1 ' 



resonances do not undergo any rescattering and thus reflect the statistics of the source resonances, which can in 
principle decrease the contribution to the HBT effect at small relative momenta. 



• S+Pb 200 GeV/A (NA44) 

long lived resonances 
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where q is the momentum transfer in the scattering, and the substitution Z — > Z + 1 takes into account 
scattering by atomic electrons p8| . While the total cross section is a — AitZ{Z + l)a 2 a 2 , the effects of 
multiple scattering are more accurately controlled by the transport cross section, at = J dfl(l — cos9)da/dfl, 
given here by 

4irZ(Z + l)a 2 „ ,„ , , 
a t = { J 1 In 2pa - 1/2 . 65 

In a single scattering, (cos 9) = 1 — at /a, so that at high energies, (9 2 ) ~ 2at/a. 

Consider a particle going a distance L through a medium of atomic density n a . From the multiple 
scattering equation one readily finds that the particle emerges with (cos#) = e~" aL<Tt , and more generally, 
with (P £ (cos6»)) = e ~ n " I " Tl , where a t = J <2f2(l - P £ (cos 9))da/dfl (39). The underlying angular distribution 
is more complicated, but for our present estimates we may assume that the spread in angles is Gaussian: 
f(0,L)9d9 ~ e~ e 'I^ 2 ^9d9. The mean square scattering angle (9 2 ) is given by the non-linear equation fl38[ , 

4nZ(Z + l)a\ aLln{e2) = ^ ^ (6g) 



P 2 



-n 



where 



a 



cff _ 4ttZ{Z + l)a 

p2 



2 



\n(pa/is 1/2 ), (67) 



and the factor v w 1.32(1 + 3.34Z 2 a 2 ) includes corrections to the cross section beyond the Born approx- 
imation. The second term on the left in Eq. (^6|) approximately takes into account effects of large angle 
scatterings, and reduces (9 2 ) from (9 2 )o. 

Let us consider as illustration the effect on a 4 GeV pion scattering through 1 mm of Pb. Then the mean 
scattering angle is ~ 2 x 10 -3 , which produces a mean transverse spread in momentum, Ap^, of order 8 
MeV, corresponding to a transverse deflection, Ar_L, of order 3 cm. Similarly, a 4 GeV pion traversing 15 m 
of air [Z = 7) undergoes a mean angular deflection ~ 0.9 x 10 -3 , with Ap± ~ 3.7 MeV and Ar± ~ 1.4 cm. 
The pion makes some 10 3 scatterings per cm; air is remarkably opaque to pions. 

The effect of these small angle deviations is to spread out the singles distributions and the correlation 
function. Such secondary scattering effects arc generally accounted for in the estimated experimental mo- 
mentum resolution. An initial distribution of single particle transverse momenta, n% ± , will be spread into a 
final distribution, 

Similarly, two particles starting out with a given relative momentum and undergoing random walks do not 
end up with the same final relative momentum. For example, two particles detected with zero relative 
momentum may have actually started out at a larger momenta and have been bent in by the air or target. 
An initial HBT correlation function Co(q) = 1 + R ° will be spread into an observed correlation function, 

C ohs {q) = l + A off e-« 2 ^«, (69) 

where the measured radius i? c ff is decreased from the original radius, R, by a factor 

R cS 1 



R [^ + 2(RA P± ) 2 } 1 / 2, 
and the chaoticity parameter is reduced from unity to 



(70) 
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For example, with an initial nominal radius of R = 7 fm, scattering in air reduces the measured R by 2% and 
the measured A by 3%. Including the effect of a 1 mm Pb target, one finds a 7% reduction in the measured 
R and a total reduction of A due to secondary scattering of 16%. The effects on HBT of secondary scattering 
in a thick target can be substantial; e.g., for 1 cm of Pb, A falls below 0.4. 

The astute reader may at this point have noticed a contradiction between the picture of secondary scat- 
tering in a cloud obscuring a smaller source, the sun say, leading to a narrower correlation function than 
the one that would be produced by the sun (cf. Fig. 7b), and the present picture of scattering in air, which 
broadens the correlation function and reduces it at the origin. I leave the resolution of this problem as an 
instructive exercise. 
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Fig. 13. E877 data @ for (a) 7T+7T , (b) 7r p, and (c) n + p, as well as comparison (dotted lines) with the toy 
model, Eq. (|78|) for tq in the range 3-15 fm, assuming a bare correlation function Co = 1, and the Gamow 
correction (solid line). 

IX. FINAL STATE COULOMB INTERACTIONS 

Up to now we have discussed Hanbury Brown-Twiss interferometry assuming that the particles travel 
completely independently once they leave the collision region. In fact one measures correlations primarily 
among charged mesons; the Coulomb interactions between any given pair of particles, as well as those of 
the individuals in the pair with the other charged particles in the system, produce important effects on the 
measured correlation functions. Even though the detectors are many meters from the collision, those that 
produce the enhanced signal are typically within a meter of each other. The pions whose correlation one 
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measures travel essentially along the same direction and they continue to have a strong Coulomb interaction 
the entire way. Disintangling these final state Coulomb interactions is a very difficult question.^] 

Bare data, uncorrected for Coulomb interactions, brings out the situation very clearly. Figure 13 shows 
the uncorrected E877 correlation function measurements for 7r + 7r~, n~p, and 7r + p pairs produced in Au+Au 
collisions at the AGS at 10.8 GeV per nucleoli jilj, while Fig. 14 shows bare data for the tt + it + and -k~it~ 
correlation functions. The correlation function for distinguishable particle pairs does not have the expected 
value of unity, while the correlation function for identical particles does not rise up anywhere as high as in 
the Coulomb-corrected data, e.g., Figs. 1, 2, and 9-11. Rather, the data for both identical particles and 
oppositely charged non-identical particles are very similar. 

The traditional method of correcting for final state Coulomb interactions is to employ the Gamow cor- 
rection. For example, in the beta decay of a neutron into a proton, electron, and antincutrino, the proton 
and electron are produced in a relative Coulomb state. Because the electron and proton are attracted 
to each other, the amplitude for them to be at the origin is enhanced. The decay amplitude is a bare 
matrix element times the relative electron-proton Coulomb wave function, -0(7(0), at the origin. Non- 
relativistically, the Coulomb wave function for the relative motion of a pair of particles of charges ze and 
z'e with relative momentum Q — (p — p')/2 = q/2 at infinity, and relative velocity v rc \ = Q/m te &, is 
V'cf^) = V , c(0)i-Fi(— irj] 1; i(Qr — Q ■ r)), where the dimensionless parameter rj(Q) is given by 

zz'a 

V = r , (72) 

Vrel/C 

and the reduced mass m rc< \ equals m/2 for two particles of mass m. Also, 

*c(0)=(^Y", (73) 



The actual rate is that which one would measure in the absence of any Coulomb effects times the Coulomb 
correction, | V'c (0) | 2 . For particles of opposite charge the probability is enhanced, by a factor tending to 
2n\r)\ at small Q. On the other hand imagine a decay of a A ++ into a proton, positron and neutrino. The 
proton and positron repel each other, and thus have a reduced amplitude to be at the origin. The Coulomb 
correction, | V^c (0) | 2 , is less than unity in this case, and the net rate would be suppressed from its value in 
the absence of Coulomb interactions, by a factor tending to 2-Kr\e~ 2 ' Kri at small Q. 

In making a Gamow correction in heavy-ion collisions, one assumes that the pair of identical particles is 
produced in a relative Coulomb state at zero separation, and thus the amplitude for doing so is reduced 
from the bare amplitude by the factor %j)c (0). The bare correlation function, Co(q), where q — 2Q, is thus 
extracted from the measured correlation function, C(q), by dividing out the assumed factor |?/>c(0)| 2 in the 
production rate: 

C(a) ( e 2*v(Q) _ i) 

Note that as the relative momentum goes to zero, r\ — » +oo and ipc(0) for identical particles, and Eq. 
((74|) yields an infinite correction. 



8 Most of the material described in this section was developed by the author and P. Braun-Munzinger, with J. Popp's 
helpful assistance, and reported in Ref. |40| . 
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The question is why one should assume that the Coulomb wave function at the origin should control the 
Coulomb corrections? For particles of the same charge, ipc(f) falls to zero exponentially as the particles 
approach each other inside the (zero angular momentum) classical turning point, r t , defined by g 2 /2m rcc i = 
e 2 /r t . Outside r t it oscillates, and describes essentially classical physics. In order for the physics at the origin 
to be relevant, it is necessary that the source be highly localized compared to the distance to the turning 
point. However, for pions in a heavy-ion collision, r t ~ (200 fm)/Q 2 , where Q is measured in MeV/c; for 
Q 10 MeV/c, a typical minimum value, r t is only 2 fm, and smaller for larger Q. Since r t is much smaller 
than the characteristic heavy-ion radius, most of the pairs of particles observed in a heavy-ion collision are 
in fact made at relative separations well outside their classical turning points. 
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Fig. 14. E877 data 41 for tt + tt + and (b) ir n , together with the Coulomb correction, Eq. ([j 
same range of r as in Fig. 13. 



for the 



There are three relevant length scales in the Coulomb problem, the classical turning point, r t , the wave- 
length of the relative motion, and the two-particle Bohr radius, ao — l/rn ro de 2 (which is 387 fm for irir and 
222 fm for 7rp). For typical Q, these length scales are cleanly separated: 



r t : l/Q : a = 2 : a Q : (a Q) 2 . 



(75) 



For 7T7T (or 7rp), a$Q — l/\rj\ — 1.96 (or 1.13) Q/(Mev/c) ^S> 1. The classical turning point is thus the 
relevant scale for Coulomb effects. These arguments suggest that the Coulomb corrections are dominated 
by classical physics. 

The major effect of the Coulomb interaction between the particles in the pair, at distances large compared 
with r t , is to accelerate them relative to each other. Particles of the same charge are accelerated to larger 
relative momenta, thus depressing the observed distribution at small Q, while particles of opposite charge 
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are reduced in relative momentum in the final state, which builds up the distribution at small Q. Although 
these effects are qualitatively similar those produced by the Gamow correction, they are quantitatively rather 
different. 

In the presence of many produced particles, the relative motion of the particles in the pair is strongly 
affected by their interactions with the plasma of other particles. The mutual Coulomb interaction of the 
pair becomes dominant only when the pair has sufficiently separated from the other particles in the system 
that there is small probability of finding other particles between the particles in the pair. 

One can write down a simple toy model to take these effects into account, by simply neglecting the 
Coulomb interaction between the pair for separations less than an initial radius ro, and for separations 
greater than ro including only the relative Coulomb interaction. Since the relative motion is in the classical 
region, conservation of energy of the pair implies that the final observed relative momentum Q is related to 
the initial momentum of the pair Qq at ro by 

£—£-±1 (76) 
2m rod 2m rod r 

where the upper sign is for particles of like charge, and the lower for particles of opposite charge. For 
example, for pions with ro = 10 fm, Q 2 = Qq ± 20(MeV/c) 2 . The physics can be treated non-relativistically 
because in the rest frame of the pair, one is interested in relative momenta small compared with mc. 

Since the Coulomb interaction conserves particles and the total momentum of the pair, the final distribution 
d 6 n/d 3 pd 3 p' of relative momenta Q is thus given in terms of the initial distribution of pairs, d 6 n° /d 3 pod 3 p' , 
by 

d 3 pd 3 p' ^ d 3 pQd 3 p' ^° ^ ^ 

The Jacobian, with changes in relative angles ignored, is, from Eq. (f76[), d 3 Qo/d 3 Q = Qo/Q. Neglecting to 
good accuracy the effects of Coulomb interactions on the singles distributions we have 

C(q) = ^C (q )= fl T ^^y /2 C7o(<fo). (78) 



Figure 13 compares the predictions of the toy model, Eq. (78), with the E877 data for 7r + 7r~, ir~p, and 
ir + p systems in Au+Au collisions at the AGS jy], assuming that the bare correlation function Co equals 
unity. The dashed lines are the results of the toy model for r = 3 fm (rightmost curve), 9 fm, and 15 
fm (leftmost curve), along with the standard Gamow correction (solid line). Except at very small relative 
momenta Q^SlO MeV/c, where effects due to the finite momentum resolution of the experiment become 
visible in the data, the model gives a good account of the data for ro in the range of 9 - 15 fm. By contrast, 
the Gamow factor considerably overpredicts the data for all Q shown here. Similar bare data from NA49 
at CERN, for 7r + 7r~ produced in 160 GeV per nucleon Pb on Pb collisions is also well fit by the toy model 
with ro ~ 10 — 20 fm, while again the Gamow correction is too large, as in Fig. 13 [42]] . 

Note that the raw correlation data for non-identical particles contains information about the mean sepa- 
ration of pairs when screening effects become negligible, summarized in the toy model by the parameter ro, 
which is possibly Q dependent. 

With the initial radius r extracted from the unlike-sign data, one can then construct the Coulomb 
correction for like-sign particles. The Coulomb correction factor deduced from Eq. (|7^) is shown for tt + tt + 
in Fig. 14a and 7r~7r~ in Fig. 14b, for the same range of ro (as in Fig. 13, the rightmost curve corresponds 
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to r = 3 fm). Again we see that use of the Gamow factor implies a correction which differs significantly 
from that of the toy model. 
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Fig. 15. Toy model calculation of C(Q) for like-sign pions (crosses), compared with the correlation function 
derived by making the standard Gamow correction (vertical bars); (a) tt + it + and (b) tt~tt~ . 

Dividing the raw E877 data by the toy model correction factor, with r = 15 fm, we obtain the correlation 
function for like-sign pions (crosses) shown in Fig. 15a for tt + tt + and Fig. 15b for ir~ir~, which also 
shows the correlation function (vertical bars) derived by making the standard Gamow correction. Using the 
Gamow factor instead of the proper Coulomb correction leads to a correlation function which is ~ 30% wider, 
implying a correspondingly reduced radius parameter. Furthermore, the shape of the "Gamow-corrected" 
correlation functions have considerable non-Gaussian tails in the range 30 < Q < 80 McV/c. These tails do 
not exist in the raw correlation function and obscure the interpretation of the data. 

The length rg gives one a measure of the scale at which Coulomb interactions between the particles in the 
pair dominate their relative motion. To calculate this decoupling length from microscopic models requires a 
non-trivial description of many-particle screening in the high frequency reqime. The Coulomb corrections will 
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furthermore change character at RHIC energies, where the meson density in given events will be sufficiently 
large that the Coulomb corrections remain those of a many-particle system out to much larger distances. 

In order to make a bridge between the toy model and the Gamow correction, as a first step in constructing 
a more accurate accounting of Coulomb corrections, it is instructive to review how the classical description 
emerges from the full quantum-mechanical treatment of the Coulomb problem. In the absence of Coulomb 
interactions (denoted by here) the number of pairs of relative momentum Q is given by Eq. (p4) : 



i/ dxdx'dx"dx'"e lP ( x+x ''- x '- x '''V 2 e^ x - x ''- x '+ x '''\j^{x)J^^ (79) 



To take into account the Coulomb interaction only between the pair of produced particles, we simply replace 
the relative free-particle wave functions, e l Q( x ~ x ) and e l Q( x ~ x ) , by the Coulomb wave functions, ipc(%— x ") 
and ipc(x' ~ x "') for the relative motion for a pair of relative momentum Q at infinity, so that 



p p d 3 P d 3 P ' 

= \j dxdx'dx"dx"'e^ x + x "- x '- x "'^ 2 ^ c (x - x")ip* c (x' - x"')(.P(x)J^(x")J{x"')J(x')). (80) 

Pairs of low relative momentum have relatively low angular momentum, e.g., a pair produced at 10 fm 
separation with relative momentum 20 MeV/c can have at most one unit of relative angular momentum. 
Thus only the low partial wave components of the Coulomb wave function enter Eq. (80) with appreciable 
probability. Let us consider just s-waves in the WKB approximation, which is quite good for the s-wave 
outside the collision volume. [] Outside the classical turning point the s-wave is 

^(r)^ rjb(r)1 ; aql ^ dn^(r), (81) 

where the local relative momentum, measuring the rate of change of phase, </>, of the wave function, is given 
by 

(Equation (|Sl|), with i'-dependent 4>(r) holds as well for higher partial waves, I > 0.) The normalization of 
( |Sl| ) agrees with ([73]) as r — > 0, while as r — > oo, the Coulomb wave function behaves as 

1>c(r) jr- sin(Qr- r?ln 2Qr + S ). (83) 

The distribution of s-wave pairs in the absence of Coulomb interactions is 

d 6 n \ s 1 



"(^) = i/ <iIdIVI " <,I '" e ' PlI+I """ ,,2 



Q\r-r"\ Q\r"-r"' 



(84) 



9 The condition for validity of the approximation is \dp(r)/dr\ p(r) 2 , which for r <C ao, the region of interest, 
becomes the restriction, r^,3/g 3 ^ 2 aj //2 . For tttt (or 7rp) pairs with Q > 20 MeV/c, WKB is reasonable for r down to 
~ 5 fm (or ~ 6 fm). 
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Then since in the region of any radius r outside the turning point the Coulomb wave function behaves locally 
as a free particle s-wave of momentum k(r), the s-wave pair distribution function is given by 



\ d 3 pd 3 p' J Q \ d 3 pd 3 p' J 



(85) 



where the factor k(ro)/Q arises from the denominators in Eq. (|8|) and (|84|). Consequently, C(q) ~ 
Co(k(r ))k(r )/Q, the result in Eq. ( |78l ) with Q = k(r ). Doing classical physics using Coulomb wave 
functions is again using a steam roller to crack a nut. 

With the connection between the toy model and the Coulomb wave function we can now extend the 
description of Coulomb corrections to smaller values of the source radius ro. In general, the effect of the 
Coulomb interactions depends on the detailed structure of the source correlation function; let us describe 
the localization of the source correlation function (,P (x)J^ (x")J (x 1 " )J (x 1 )) in both \r — f" | and |r ' — r'" \ 
by writing 

{J\x)JHx")J(x"')J(x')) m s(x - x")s(x' - x" 1 ) 
x «jt(z) J(x')) (ji(x")J(x'")) + (j\x)J(x"')) (j\x")J(x')} ), (86) 

where s(x— x') defines the effective initial separation of the pair. For small relative momenta, the combination 
ipc{x)s{x) = f(x) varies slowly on a scale of the coherence length £ c . Then we find, roughly, 



n 



£p£p d 3 pd 3 pi 

= U dxdx'dx"dx'"f{w){f*{w) + f*(-w)y p( - x+x ''- x '- x '''V 2 (jl(x)J(x'))o(Jl(x")J(x'")}o, (87) 

where w — (x — x 1 + x'" — x")/2. The integrals in this equation are sufficiently involved that the Coulomb 
corrections would have to be extracted numerically. However, if as an approximation we simply replace the 
term f(w)(f(w) + f(—w)) by its integral over all space, we arrive basically at Pratt's formula E3]: 



C(Q)~ J d 3 r(\f(r)\ 2 + f(r)r(-r))C Q (Q), 

for the modification of the correlation function by Coulomb interactions. 
The correction to the direct term in Eq. (J88|) has the form 



C(Q) dir = J IVcMIXr)!'. (89) 

To illustrate the transition from the Gamow correction to the toy model let us take |s(r)| 2 to be a normalized 
Gaussian of range ro'- \s(r)\ 2 = (27r)~ 3 / 2 r^ 3 exp(— r 2 /2r^). We show, in Fig. 16, for the -k + tt~ system, the 
results of calculations of the correction term J \ipc{i")\ 2 l s ( r )| 2 f° r r o — 1; 5, 9, and 18 fm (dash-dot curves, 
the highest for ro — 1 fm, and falling with increasing ro). As ro — * 0, the projection of the square of the 
Coulomb wave function onto the source |s(r)| 2 converges to the standard Gamow correction (solid line); for 
ro < 0.1 fm (not shown in Fig. 16) the difference between the Gamow correction and a calculation with 



Eq. (89) is less than 0.5%. For larger ro values the correction rather quickly approaches the prediction of the 
toy model (shown here for an initial radius of 9 fm as a dotted curve), indicating that, for pairs originating 
outside their classical turning point, the toy model provides an adequate and reasonably accurate description 
of the Coulomb effects. 
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COULOMB CORRECTION 
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Fig. 16. Transition from the toy model (dotted line, with ro = 9 fm) to the Gamow correction (solid line) 
with decreasing source size, calculated from Eq. ( |S9| ) (dash-dot curves). From highest to lowest dash-dot 
curves the source range ro is 1, 5, 9, and 18 fm. 

Let us turn next to the question of the effects of the Coulomb interactions of the pair with the remaining 
particles. This is a difficult many-body problem, which we greatly simplify as a first approximation by 
assuming that the remaining particles can be described by a central Coulomb potential, Z c ffe 2 /r, where in 
a central collison of nucleus A with nucleus B the effective charge Z c s is of order of the total initial nuclear 
charge (Za + Zb)- This central potential accelerates positive mesons away and slows down the negatives, 
effects described by the Coulomb wave functions for the potential. The final momentum of any particle is 
related to the initial momentum p a at production point r a by 

e (p)=e(p a )±^-, (90) 

Ta 

where e(p) = (p 2 + m 2 ) 1 / 2 . (While Coulomb effects for the relative momentum can be treated non- 
relativistically as in Eq. (|76|), the individual momenta are generally relativistic.) For simplicity let us 
ignore quantum mechanical suppressions or enhancements of the amplitude for particle emission, as well as 
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possible effects of angular changes in the individual particle orbits on the particle distributions. Then the 
single particle distribution is modified by the central potential, analogously to Eq. (fn\), by 

d 3 n(p) _ d 3 n (p a ) d 3 p a _ p a e{p a ) d 3 n (p a ) 
dp 3 d 3 p a d 3 p pe(p) d 3 p a 

Both the magnitude of the distribution as well as its argument are shifted. Experimental observation of 
these effects is reported in Ref. Q. 

Although the central potential shifts the singles distribution, it cannot introduce any correlations among 
emitted particles that have no initial correlation in the absence of the central potential, e.g., as one usu- 
ally assumes for different species or oppositely charged pions. If in the absence of the central potential, 
uncorrelated particles [C{Q) = 1] arc emitted in independent free particle states, then in the presence of 
the potential they are emitted in Coulomb states for the central potential, but still d 6 n(p , p') / d 3 pd 3 p' = 
(d 3 n(p) I d 3 p)(d 3 n(p') / d 3 p') and C(Q) remains unity. 

For particles that are initially correlated as a consequence of Bose-Einstein statistics, d 6 n(p , p') / dp 3 dp' 3 
and (d 3 n(p) /dp 3 )(d 3 n(p') /dp' 3 ) will be modified both by the Jacobians of the transformations from initial 
to final momenta, and shifts of argument. However, in forming C(q), the effects of the Jacobians in the 
numerator and denominator essentially cancel, and the primary effect is the shift in the arguments: 

rr , _ {d e n 2 (p a ,p a ')d 3 Pa d 3 p' a } 

° WJ {(d 3 n(p a )/d 3 p a )(d 3 n(p a ')/d 3 p' a )}- iy > 

Since positive particles are accelerated, the final momentum difference, q = p — p' , of a positive pair will 
generally be larger in magnitude than it is initially, while for negative pairs the final momentum difference 
will generally be smaller. Thus we expect the central Coulomb potential to cause the size of the collision 
volume extracted from positive pairs to be smaller than the actual size, and that from negative pairs to be 
larger than the actual size. As an illustration consider a pair of rclativistic particles whose initial momenta 
p a and p a ' are equal in magnitude to p a , and final momenta p and p' equal in magnitude to p\ then 

q=(p/p a )q a ^q a (l± Zc ^l ra y (93) 

where the upper sign refers to both particles positively charged and the lower to both negatively charged. 
For Z ~ 150, r a ~ 7 fm and p a ~ 300 MeV/c, the effect is an increase for positives (and a decrease for 
negatives) in the observed scale of C(Q) and decrease (or increase) in the extracted radius of ten percent. 
Such a shift of the same magnitude has been observed by E877 in 10.8 GeV/A collisions of Au on Au [ [TT| ; 
however, NA44 recently reports an effect in the opposite direction, in 158 Gev/A Pb on Pb collisions, in 
radii as a function of charged particle multiplicity |?3| , indicating the need for a more refined theory of the 
effect of the central Coulomb potential pa]. 



X. APPLICATIONS IN CONDENSED MATTER AND ATOMIC PHYSICS 



Let me finally mention briefly work on HBT in condensed matter and atomic physics. Recently, Yasuda 
and Shimizu (at Tokyo University) have made the first measurement of HBT correlations in an atomic 
system, observing the time correlations in laser-cooled ultracold (but not yet Bose-Einstein condensed) 
beams of bosonic 20 Ne atoms. The correlations in the beam are those expected from a thermal source, where 
the correlation time is the inverse of the temperature of the beam. Indeed the HBT correlation function 
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begins to rise at time separations less than ~ 0.5 x 10~ 6 sec to a value a factor of two larger than at large 
time, corresponding to a beam temperature ~ 10 2 /iK. This result is very similar to the original Hanbury 
Brown-Twiss tabletop experiment on photon bunching from a Hg vapor lamp. By contrast, a measurement 
of HBT correlations in the MIT atomic laser |}5| would yield no such atomic bunching, because of the 
coherence of the beam, but rather the correlation function would remain flat. Lack of an HBT enhancement 
would indicate coherence of the beam. In general, loss of HBT correlations would probe the onset of Bose- 
Einstein condensation, not only in atomic systems, but in condensed matter systems such as the observed 
Bose-condensed paraexcitons in cuprous oxide jl6|,fl7j . 

Another interesting application of HBT has been in light scattering from atoms trapped in optical lattices 
. Jurczak et al. (at Orsay) have created an optical lattice with an arrangement of four lasers in which 
they trap atomic rubidium at a density ~ 2x 10 9 cm~ 3 , filling about 10~ 4 of the lattice sites. The lasers also 
scatter from the rubidium, and the time correlations in the scattered light (of two different polarizations) 
effectively measure the atom-atom correlation functions in the lattice. From these measurements they are 
able to measure the diffusion of the loosely packed atoms in the optical lattice. Lastly we mention that HBT 
has also been proposed as a probe of the space and time structure of bubbles in sonoluminesence Q . 

In summary, the technique of Hanbury Brown and Twiss, which was first developed to measure astronom- 
ical object of sizes at least 10 12 cm, has, as we have seen, turned into a valuable tool to measure subatomic 
phenomena on the quite opposite scale of 10~ 12 cm. More recent experiments have shown its utility in 
atomic and condensed matter physics as well. While the basic theory underlying the nuclear applications is 
established, as described in these lectures, many effects, e.g., Coulomb interactions, possible non-chaoticity, 
non-zero coherence lengths, multiple scattering, etc., introduce various levels of uncertainity into the in- 
terpretation of the HBT measurements. A better understanding of such effects remains a challenge in an 
accurate connection of HBT measurements to the microscopic physics of collisions. 

These lectures are a small birthday tribute to my dear friend Wieslaw Czyz who over the years opened 
many worlds to me - from Zakopane and Cracow, to the pleasures of high energy nuclear physics. I would 
like to take this opportunity to express my gratitude to the organizers of the present Cracow School of 
Theoretical Physics at Zakopane for enabling me to participate in the School where these lectures were 
given. I would also like thank the members of my group in Urbana - Alejandro Ayala James Popp, and 
Benoit Vanderheyden - who are responsible for much of the material reported here, and Michael Baym for 
preparing the graphics. I am also grateful to Peter Braun-Munzinger, Henning Heiselberg, Barbara Jacak, 
and Dariusz Miskowiec for many discussions of this material and for making figures available, and to Ulrich 
Heinz for insightful comments on the manuscript. This work was supported in part by U.S. National Science 
Foundation Grant No. PHY94-21309. 
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